TAUP 2668-2001 



Relativistic Mechanics of Continuous Media 

S. Sklarz 

Department of Chemical Physics, Weizmann Institute of Science, 

Rehovot 76100, Israel 
shlomo . sklarz® weizmann. ac .il 

L. P. Horwitz 
Raymond and Beverly Sackler Faculty of Exact Science, 
School of physics, Tel Aviv University, 
Ramat Aviv 69978, Israel [] 

In this work we study the relativistic mechanics of continuous media on 
a fundamental level using a manifestly covariant proper time procedure. We 
formulate equations of motion and continuity (and constitutive equations) 
that are the starting point for any calculations regarding continuous media. 
In the force free limit, the standard relativistic equations are regained, so 
that these equations can be regarded as a generalization of the standard pro- 
cedure. In the case of an inviscid fluid we derive an analogue of the Bernoulli 

equation. For irrotational flow we prove that the velocity field can be derived 
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from a potential. If in addition, the fluid is incompressible, the potential must 
obey the d'Alembert equation, and thus the problem is reduced to solving 
the d'Alembert equation with specific boundary conditions (in both space 
and time). The solutions indicate the existence of light velocity sound waves 
in an incompressible fluid (a result known in previous literature 0]). Re- 
laxing the constraints and allowing the fluid to become linearly compressible 
one can derive a wave equation, from which the sound velocity can again be 
computed. For a stationary background flow, it has been demonstrated that 
the sound velocity attains its correct values for the incompressible and non- 
relativistic limits. Finally viscosity is introduced, bulk and shear viscosity 
constants are defined, and we formulate equations for the motion of a viscous 
fluid. 

1 Introduction 

We wish to formulate manifestly covariant equations for the motion of a 
continuum on a fundamental level, and to demonstrate their use in an inter- 
esting problem. In attempting to reach a consistent theory one meets a few 
difficulties: 

The general motion of a continuum may include accelerations, making 
the use of standard special relativistic dynamics problematic. The usual pro- 
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cedure for dealing with an accelerating particle, would be to go from one 
inertial frame to the next, as the particle accelerates, always keeping it in- 
stantaneously in an inertial frame In the case of a continuum however, 
this procedure is impossible, for each of the infinitesimal volume elements ac- 
celerates, in general, at different rates and directions, and so only one specific 
portion of the volume can be kept at rest at any particular time. Secondly 
the above procedure may not be manifestly covariant, as an accelerating 
particle momentarily at rest in some frame, is obviously not equivalent to 
an identical particle moving at constant speed. The accelerating particle, 
for instance, will be radiating if it is charged. Our claim is that keeping a 
particle momentarily at rest can, at best, yield an approximation to the real 
motion, but is by no means a consistent covariant treatment of the system. 

Another severe difficulty in writing down relativistic equations of hydro- 
dynamics lies in the relativistic property of simultaneity. In order to define 
an extended body in mathematical terms, a certain configuration of the body 
is chosen arbitrarily and is referred to as the reference configuration. In the 
nonrelativistic theory, each particle, or infinitesimal volume element of the 
body, can then be labeled throughout the evolution by its position in that 
reference configuration. The difficulty in constructing a relativistic theory 
of hydrodynamics lies in the fact that each inertial system will have its own 
notion of time and simultaneity, and they will not be able to agree on any set 
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of particles to be considered as existing simultaneously at any given time t. It 
therefore becomes impossible to define a consistent reference configuration. 

In 1941 Stueckelberg [|| postulated that events (not particles) are to 
be considered as the fundamental dynamical objects of the theory of mo- 
tion. The dynamics of an event is governed by 8 independent parameters 
consisting of the space-time coordinate^ = {ct, r) and energy momenta 
pf^ = [E/c, p) associated with each event, and parameterized by an invariant 
time. Based on this, Horwitz and Piron in 1973 proposed that there exists 
a universal time r, by means of which dynamical interactions are correlated 
||13|| . This universal time may be identified with the Robert son- Walker time 



of the expansion of the universe, the time read on a freely falling ideal clock. 
The evolution of a system of events through r is generated by an invariant 
function K on this 8A^ dimensional phase space, by means of the covariant 



Hamilton equations ||11|| : 



dxf OK dpt OK 

— = ) = ) (1) 

dr dpi^' dr dxi^' 

The collection of events along each world line corresponds to a particle, and 
hence the evolution of the state of the A^-event system describes, a posteriori, 



^Here and throughout this work we use the usual convention where Lattin indices 
k, ... take the values 1, 2, 3, related to Euclidean space, whereas Greek indices fj,, v, ... 
take the values 0, 1, 2, 3 and correspond to the four dimensions of Minkowski space-time 
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the history in space and time of an A'"-particle system. In the case of a free 
particle, for instance, we choose for the Hamiltonian: 



where M is a positive parameter, which is a given intrinsic property of an 
event, having the dimension of mass and emerging as the Gahlean target 
mass in the nonrelativistic hmit c — > oo. This results in 



so that we have for the physical velocity (we take c = 1 in what follows unless 
otherwise specified), 

in total agreement with the usual Einstein kinematics. As {E, p) are each 
taken to be independent variables of the motion, the particles are not nec- 
essarily on mass shell, and the particle mass squared: rri^ — E"^ —p^, is not 
contrained to be constant (except for a free particle), but is rather to be 
determined as a solution of the dynamical equations. Note that (we use the 
metric — h +- 1-) 
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~~d^~d^^~d^^Jp 
and is unity only for the "mass-shell" value = M^, Hence r corresponds 
to the proper time of an ideal free clock on its mass shell. On mass shell, the 
time obeys the relation: 



dt E 1 

(6) 



dr M y/l 



V 



2 ■ 



SO that dt is precisely the time interval measured in the laboratory between 
two signals emitted by a source, traveling with velocity v, with interval dr, 
according to the Lorenz transformation. If the emitter is not on shell there is 
a factor m/M. The red shift imposed by general relativity. At = y^^Ar is 
another example of the effect of dynamics, in this case, the acceleration due 
to gravity (contained in goo), on the observed time interval. The observed 
time interval is therefore influenced by forces that may move the energy 
momentum off shell. 

This theory has since been applied to fundamental questions such as the 
Newton- Wigner position operator and the Landau- Peierls relation 0] , and to 
dynamical problems such as two body bound states ^ and scattering and 
electromagnetic interactions in classical mechanics 0, §|, quantum mechanics, 
and quantum field theory H, [10 . 
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In order to formulate equations of motion for continuous media, we there- 
fore turn to the above mentioned procedure of Stueckelberg, Horwitz and 
Piron , and all following discussions will be conducted in the framework of 
this theory. As a consequence of the existence of a universal time which pa- 
rameterizes all systems, it becomes possible to define a certain configuration 
of events in space-time X^, at some definite initial universal time tq = 0, 
and then follow the evolution of the system through r. The universality of r 
eliminates any ambiguity as to the identification of events and their distribu- 
tion in space-time. Labeling the various events by their initial configuration 
X^, one can thus express their respective positions at any later time r: 

x'^ = x^(X",r). (7) 
The velocity and acceleration of an event labeled at time r is thus: 



V.(A-.r) = ^^^(^ (8) 

CLT 

^^^'^^ = Tr = dr^ 

The Eulerian velocity field describes the velocities of events passing a certain 
point in space-time, as a function of the universal time: 



v^'{x'', r) = ViX'^ix'', r), r) (10) 
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We now wish to express the acceleration in terms of the Eulerian velocity 
field: 



A^{X'',t) 



(11) 





A''{X''{x\t),t) 



(12) 



dr 



+ 



dr dx' 



dr 



dx' 



In space-time, due to mass-energy equivalence, there is no a priori con- 
servation of mass of individual particles; rather we shall be dealing with a 
"conservation of events" . Let n{x'^, r), be the density of events (this quantity 
is a scalar since the four volume element is invariant) at the space-time point 
x^^ at time r multiplied by M, the intrinsic, given mass dimension property 
attributed to each event. Consider now a small volume of fluid. The flux of 
events through a closed surface enclosing the volume must cause a decrease 
in the density of events within the volume. It is in place to stress at this point 
that the surface mentioned is a three dimensional surface, the direction of 
which is defined by a four vector normal to that surface in four dimensional 
space-time. Hence, for example, we define dS^ as a volume of three dimen- 
sional space, its direction being the time direction. Likewise, the volume over 
which the integration is carried out is obviously the four dimensional space- 
time. The fiux of events past an infinitesimal surface dS^ is: e{v)nv^dS^, 
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where e{v) = ±1 depending as to whether is space-hke (+1) or time-hke 
(— 1). The reason for introducing this additional sign is that for and dS'^ 
both pointing out of the bounded volume, in a time-hke direction, the vector 



product v^dS^ is negative whereas the flux of events must be positive 119 
Therefore: 



d_ 



ndV = / eMnv^'dS. (13) 
V Js 



or, using the four-dimensional Gauss law, ||18 



^^dV = J^e{v)nv,dS^, (14) 



we have. 



,^dV^[ ^dV (15) 

Iv or Jv ax'^ 

and so taking the volume to be infinitely small: 

which is the equation of continuity for the flow field. Note that in the case 
where the density of events n is constant, the above equation reduces to: 



dX = 0. (17) 
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In a continuum we focus on an infinitesimal volume element and describe 
the forces acting on it. In general one talks about two kinds of forces: 

1. Forces which are proportional to the density of mass, charge or number 
of events, within the volume. The force will therefore be nb^, where n 
is the density of events. 

2. Surface contact forces, which are exerted by the continuum of events 
surrounding the volume under consideration, and are transferred through 
its surface. The number of events neighboring a small portion of the 
surface, dS, in a locally homogenous fluid, are obviously proportional 
to its size, hence we can assume the force to be proportional to the sur- 
face it is acting upon. In mathematical terms we assert that the force 
df on an element of surface dS tends to zero as dS tends to zero, but 
that the ratio df /dS tends to a definite limit. This assumption enables 
us to define a stress four tensor, the components of which describe the 
force in the n direction, exerted on a surface in the u direction: 



The surface force on an elementary volume unit, in the ii direction is 
therefore (by the usual argument): 

nr = ^ (19) 
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Summing up the forces acting on an infinitesimal volume element and 
equating to the density of acceleration within the volume, we get the equa- 
tions of motion: 

By introducing a manifestly covariant Boltzmann equation, Horwitz, Shashoua 



and Schieve [0 were able to write down equations for conserved quantities 
in a statistical mechanical framework. It is interesting to note that the equa- 
tions of continuity and motion derived above from a hydrodynamical point of 
view are equivalent to those presented by these authors for the conservation of 
some constant quantity and for the conservation of momentum respectively. 

In order to make these equations of motion operational it is necessary to 
formulate constitutive equations for the fluid. These will define the way by 
which the stress forces acting at a certain point in the fluid a'^" , are governed 
by its physical state (and/or history). These laws however are determined 
by physical properties of the material; we shall make some simplifying as- 
sumptions about the nature of the fluid, comparable to those which lead to 
the Navier-Stokes equations in the non-relativistic case. 

In the following work we shall formulate a potential theory for inviscid, 
irrotational, incompressible fluids. Then, relaxing the constraints, we shall 
allow the fluid to become linearly compressible, an assumption which will 
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lead to acoustic modes within the flow. Finally we shall introduce equations 
for the motion of viscous flow. 



2 Standard relativistic limit 

In the current section we shall show how Eq. ([16|) and (0) transform into 
the standard relativistic equations of hydrodynamics, which in turn give the 
nonrelativistic limit by taking c — »• oo. We begin with a few preliminary 
remarks which will assist us in what follows. Note the following connection 
between the relativistic force term and its nonrelativistic counterpart P: 



M 



M 



M 



dx^ d 



I dx^' 



' dx^ dx^ ' 
dx'^ dr 



dx^dh^ 
dx^ dr"^ 

dx^ (P'x^ 



dr dr \ dx^ 



+ 



'd^ 
dr 



+ 



(21) 



Taking the space components we can write 



r = 1-1 f' + v7 



:22) 



A rearrangement of terms gives the physically measureable force f* in terms 



of the relativistic force, 



12 



f^= (^)'(f -v7°). (23) 

Another useful expression which connects derivatives of the velocity field 
with its nonrelativistic physically observable counterpart v* = cfYf", will be 
used later and can be proved as follows (we denote by ds a derivative with 
respect to any variable): 

d,v' = ds{ ) = -dy + (24) 

c c c 

We remind ourselves too, that v*^ = dx^/dr = E/M, and that in the 
nonrelativistic (on-shell) limit f ° = c/'^Jl — ^ ~ c(l + |^). 

To put our equations into correspondence with standard results, we must 
note that any interaction with an apparatus or some macroscopical object 
corresponds to an interaction with a worldline, an integration over r is 
therefore necessary. In performing this procedure to obtain the on shell 
limit, we assume that correlations in r survive only over a short inter- 
val, thus a term of the form: / a{T)h{T)dT can be factored to obtain 
~ / a{T)dT J b{T')dT' = ab in a T averaged sense. 0. 

^ If both air) and 6(r) have support only arround the zero mode then 
J a{T)b{T)dT ^ J a{s)b{-s)ds 
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Applying this procedure to equations (p!6| ) and (^) causes the r derivative 
terms in both equations to drop due to the assumption that both n and f " 
vanish pointwise as r — cxd. From Eq. (^) we then get: 







do{nv^) + di{nv^) 



d_ 

di 



(26) 



V 



,0 




,0 



(27) 



where all quantities are from now to be understood as being r averaged. 

Written in the form (^), Eq. ([T6|) can be interpreted as expressing the 
conservation of energy density (or the continuity of energy flow), whereas the 
second form (^) states that the conservation of mass is to be corrected by 
a relativistic term. This can be seen by noting that the first bracketed term 
in (P7|), is just the regular nonrelativistic equation of continuity which gives 
the change in time of the mass of a small volume of three dimensional space. 
In nonrelativistic hydrodynamics this term should obviously be zero, for the 
mass is conserved. This however is not the case for a relativistic fluid, as an 
increase in the energy also has to be taken into account in the conservation 
law. This is precisely the meaning of the second term. It expresses the change 




(25) 
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of energy of the small volume of fluid moving with the particle. This change 
of energy causes a corresponding non vanishing change of mass which must 
balance it out in such a way that the total be zero. It should be noticed 
that the derivatives of f ° which appear in the second term of (|27| ) are in the 
standard relativistic limit, of the order 1/c, and so are small compared to 
the first term. Therefore when taking the nonrelativistic limit this second 
term can be neglected, and one is left with the standard equation for the 
conservation of mass: 



We shall now treat the space and time components of Eq. (pOD separately. 
We consider a fluid free of any external forces and therefore take 6 = 0. For 
the sake of simplicity of notation we assume here that a^" = —pg'^'^. A 
physical justification for this simplification will be forthcoming in the next 
section, but the generality of our following discussion is not based on or 
restricted by it. Taking the time component we have: 



dn 
'dt 



+ V ■ (nv) = 0. 



(28) 



-9° 



P 



nv 





or. 




(29) 
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This gives the time change of energy moving with the fluid. Inserting the 
standard relativistic hmit for f", and neglecting terms of order we are 
left with 



= ^(^^^') + V-(vlnt;2), (30) 

where we made use in the third line of the equation of continuity (|28|). Eq 
(|30| ) is evidently an equation for the change in nonrelativistic kinetic energy. 
The first term on the right hand side can be interpreted as the change of 
kinetic energy within a small volume of space whereas the second term is the 
flux of energy leaving that volume. 

Turning to the space components, and using relation (|2^) we get: 



= -(t;V9^v' + vV^y). (31) 
Inserting Eq. ( pQ]) and rearranging terms yields 
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By remembering relation ( p2D one can notice at once that the left hand side 
of Eq. (P^D is simply the physical force acting on an infinitesimal volume of 



fluid, whereas the right hand side gives the rate of change of the physical 
velocity moving with the fluid. In fact it can clearly be seen that in the 
nonrelativistic limit c oo, v^/c — » 1, Eq. (|3^ ) simply becomes the non- 
relativistic Euler's equation of hydrodynamics. It is worth noting too, that if 
we accept a correspondance between the density of events n and the sum of 
the internal energy and pressure w = e+p (also called the heat function), then 
Eq. (|3^) is precisely equivalent to the relativistic hydrodynamics equation 
derived by Weinberg and Landau in the forceless events framework 



of standard relativity for which (f°/c)^ = 7^. 



3 Potential Theory for Relativistic Hydrody- 
namics 

Under certain circumstances some problems in fluid dynamics can be solved 
in an approximation in which viscous forces are neglected, and the fluid is 
assumed to be incapable of sustaining shear forces. If the fluid is assumed 
also to be spatially and temporally isotropic, we can write down the stress 
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tensor as follows: 



a 



-P9 



p when fi = u = 
—p when /i = z/ 7^ 
when ^ ^ 



p 














-p 














-p 














-p 



, (33) 



where g^'^ is the metric tensor for flat Minkowski space and p is a Lorentz 
scalar. The interpretation of the stress tensor components in this form, is as 
follows: The three space-space components relate to the force exerted within 
the fluid per unit surface per unit time, whereas the time-time component 
describes the forces acting in the time direction per unit volume of the fluid 
(essentially mass-changing forces). 

Substituting Eq. (Rl) into the equations of motion. 



OT 



(34) 



this leads to 



- d^p + nb" = n(— + v''d.)v''. (35) 
or 

Eq. ( p5D may be written in different form by introducing two new variables. 
We define ^ = r:^^^dn', so that ^d^'p = 9^*. Note that for incom- 
pressible fluid n is constant and we have simply \I' = p/n. Assuming too 
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that the body force can be derived from a potential function b"^ = —d'^Q , we 
can write: 

-9^(^ + ^]) = (^ + ^;^5^K, (36) 

An important kinematic result which shall be used later, is known as the 
Circulation Theorem. Consider a closed circuit C linking a continuous line of 
fluid events. In general the space-time configuration of the loop C depends 
on r. We denote this symbolically by Cr- The kinematic Theorem to be 
proven is : 

where is the space-time acceleration field of the fluid. The proof can be 
obtained by transformation to the material coordinates as follows: 



i v,dx^ = i V,-^dX^. (3^ 
JCt J Co dX^^ 



So then: 



drJc. ^ J Co dr [dXi^J Jco dr dX'' ^ ' 

The first integral on the right hand side of Eq. (|39|) , can be shown to vanish: 
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since the path is closed. The second can be transformed back to the spatial 
coordinates: 



proving the theorem. 

Substituting for in Eq. (pT]), from the equation of motion (|36| ) gives: 



d_ 



[viz + = 0. (42) 



where it has been assumed in the last step that the functions are single 
valued. This leads generally to: 

^ Vf.dx'' = Const (43) 

If the motion of the fluid was generated from rest, so that for some initial 
time the circulation was zero, we get: 



^ v^dx'' = 0, (44) 
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for all times. A necessary and sufficient condition, for Eq. (0), with contin- 
uous f ^, is that v'^ be expressible as a gradient of a potential function. We 
write 

= -d^^. (45) 

where $ is defined as the velocity potential. It has been shown above (Eq. 
(p!7| ) that for an incompressible fluid 

ax = 0, (46) 

so combining the two, we find that the velocity potential, $, satisfies the 
d'Alembert equation: 

= 0, (47) 
providing an enormous simplification in the theory of inviscid irrotational 



flow. Using the four dimensional Stokes theorem |2T| we can convert the line 
integral in Eq. (|4^) into a spatial integral: 



v^dx^ = / / {d^v, - d,v^)dS^'' = 0. (48) 
If this is to be true for arbitrary surface S, the integrand must vanish, giving 



{d^v, - d,v^) = 0. (49) 
Using this we can evaluate the following expression: 
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v^d^v^ = v^d^v^ = \d''{v^v^). (50) 

Inserting Eq. ( pO]) into the equations of motion (^), and rearranging terms, 
we obtain: 



5-(vI/ + ^] + l(^;^^;^')) =^Vr (51) 



Inserting Eq (^), into (|5T|) and rearranging terms, leads to further simphfi- 
cation: 



_ a- ( VI/ + + ^-{vX) -^'^1=0, (52) 



which means that: 



d_ 



m + n + -{v^v^) - — $ = /(r), (53) 



where /(r) is some arbitrary function of r. This can be considered a rela- 
tivistic equivalent to the Bernoulli equations^. 



We return now to Eq. (47) and assert that under the assumptions made, 
our problem reduces to solution of the d'Alembert equation, given specific 
boundary conditions. These require that the component of the velocity field 



^By this we mean solely that this equation fulfills, the same function in the logical 
procedure as does the Bernoulli equation in the nonrelativistic counterpart. 
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normal to the boundary vanish everywhere on its surface. Thus the motion of 
the fluid will be governed solely by the space-time geometry of the problem. 
The solution can then be inserted back into the equations of motion (^3|), 
determining the distribution of forces and stresses within the fluid. 

The general solution to the d'Alembert equation (|47|), can be written as 
an integral over all possible plane waves: 

$ = y"c/^fcA(F)exp(zA;^x^), (54) 
where k'^ is restricted to the shell 

= 0, (55) 

and the A{k^) which determine the shape of the wave packet are to be de- 
termined by boundary and initial conditions. 

It must be noted though, that there also always exists another trivial 
solution, because any solution to the Laplace equation with an additional 
linear time term, will also satisfy the d'Alembert equation. The space part, 
$V) by definition, vanishes under the Laplacian, while the time component 
being linear, vanishes when differentiated twice. This trivial solution 

^NR = ^V+UoX^, (56) 

gives rise to a physical velocity field: 
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= % = -od'<l>v (57) 
Vnr w° 

which is clearly to be interpreted as the nonrelativistic flow of an incompress- 
ible fluid. The total solution then should be written as a linear combination 

$TOT = ^NR + (58) 

As for the physical motion of the fluid we now have: 

V' = 4^ (59) 
Vnr + ^° 

If we assume, as is reasonable for low energies, that the nonrelativistic part 
of the flow is far greater than the relativistic correction, v'^j^ ^ v'^, then an 
expansion in a Taylor series gives: 



_ Vnr , VnrV° 



+ -.- T^- (60) 



^0 ^0 (^0^2 



The first term is just the nonrelativistic, t independent flow of Eq. 
whereas the other two are of a smaller order and describe time dependent 
wave-like fluctuations. It seems, then, that the covariant form for the po- 
tential flow gives rise to a physical velocity which can be divided into two 
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parts V ^ ^NR + ^R, an underlying nonrelativistic velocity field and an addi- 
tional relativistic correction with a wave-like nature, due to disturbances in 
the background fiow. These disturbances can be envisioned as ripples riding 
above the usual nonrelativistic fiow, at a speed which will be shown to be 
the speed of light. 

The velocity of the disturbances can be calculated in various ways. Writ- 



ing fc'^ = kn^^ in Eq. (|5^, and picking out a specific direction n'^, we can 



perform the integration over the magnitude k, getting: 

$ = fin.x^. (61) 

But from relation (0), n^n^ = 0, so the phase velocity vj, of the wave front 
is: 



rii 



(62) 



If the solution given above (|5J), describes a wave packet, the group veloc- 
ity of this packet can be determined through the stationary phase method. 
The highest contribution to the integral will come from those fc's for which 
the phase is stationary: 
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(63) 

so, denoting by Vg, the group velocity of the wave packet, we get 













6x^ 




dk^ 



where in the last step we used relation (pSl) . 

It is in point to stress that it is not the fluid which flows with the speed of 
light, rather it is the disturbances propagating in the flow field which travel 
with the above calculated group and phase velocities. 

In previous work on the subject, there has been much debate as to the ac- 
tual definition of incompressibility of relativistic fluids. The most straight for- 
ward definition in the framework of standard relativity, would be to demand 
that an element of the fluid shall retain the same proper volume throughout 
its motion, or in other words that the expansion of an elementary world-tube 
of liquid shall be zero. The difficulty lies in the fact that this definition yields 
infinite sound wave velocities. This would violate the most basic assumptions 
of standard relativity which forbids any material particle or information to 
be propagated at a speed higher than light. A disturbance such as a sound 
wave does carry energy, momentum and information and should therefore 
not break the light speed barrier. For this precise reason Synge intro- 
duced two more alternate definitions, the second of which actually defines an 
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incompressible fluid as such that sound waves propagate within it at hght 
velocity. 

In comparison, the equations we have written down and the simple as- 
sertion that incompressibility is to be defined on events in four dimensional 
space-time, lead smoothly and naturally to light speed sound waves in agree- 
ment with Synge's criterion. Moreover it can be seen that as c ^ oo, the 
sound waves speed, tends to infinity too, as would be expected of an incomm- 
presible fluid in the nonrelativistic limit. 

4 Compressible fluid, Acoustic approximation 

We wish to widen the scope of our study to encompass the behavior of a 
compressible fluid, and we shall show that this additional freedom introduces 
acoustic modes with sound velocities other than that of light. 

We begin by assuming that the density and pressure vary only slightly 
from some constant values. 

n = no n', P = Po + p'- (65) 

Generally the pressure is a function of the density so one can expand the 
pressure in a Taylor series: 
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pin) = p(no) + n — = p(no) + u^n (66) 

on 

where higher derivatives in n have been neglected and we have in the last 
stage denoted the constant = Comparing equations ( p5| ) and ( pB| ) 

we get a linear relationship between the variations of density and pressure, 



p' = u^n'. (67) 

We wish also to assume the three space components of the velocity field to 
be, in some sense, small. The time component on the other hand must by 
definition be of the order of c and can therefore not be considered small. This 
reasoning leads us to divide the velocity field too into two parts, 

v^' = v^ + v'^", (68) 

where Vq = [vq, 0, 0, 0) is a constant^, pure, time-like vector and v'^ is a small 
perturbation the nature of which we wish to determine. It has been shown 
earlier that the motion of the fluid, neglecting unimportant body forces, is 
governed by the following two equations: 



-d''p = ni — + v''d,vn (69) 



^It is sufficient if vq varies only slowly in spacetime relative to v'^; this possibility will 
be discussed later on. 
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-nd.v^" = + v^'d.n (70) 
or 

We assume firstly, that all changes in r in equations (^) and (ffOD, are 
negligably small compared with the spacetime changes, and we can therefore 
omit the r derivatives from both equations. This can be explained too on 
grounds of integrating the above equations over r and asserting that all 
physical quantities vanish as r — ±00 or by taking the zero component of 
the frequencies in r. All three arguments are essentially equivalent and are 
based on the fact that any apparatus in the laboratory is not capable of 
resolving fast changes in r, but can rather only measure averages over large 
periods of the universal time. The resulting equations 

- d^p = nv^d^v^, (71) 
-ndf.vf' = v^'d^n, (72) 

can now be linearized by inserting relations (^) and (^) and keeping only 
terms of the first order (n', p' and v'^ being of the first order). Higher order 
terms are neglected and so in this approximation the two equations reduce 
to 

- u'^d^'n' = nov^d^^", (73) 
-nod.v'^ = v^d^n', (74) 
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where by remembering relation (p7|), we have ehminated %>' from equation 
([7T|). The importance of the background vector fg, now becomes apparent. 
It ensures that the velocity field remains time-like and that the large time-like 
component does not get neglected in the linearization process. 
By taking space derivatives of Eq. ( [73D it can be shown that 

vld„{d^v, - d,v^) = 0, (75) 

implying that if the flow was initially irrotational everywhere, then it will 
remain so throughout the entire evolution and that the velocity field can be 
derived from a potential: 

v^' = -d^"^. (76) 
Multiplying Eq. ( [73D by Vq^ one obtains 

- u^vi^d^n = noVof^Voduv'^ , (77) 

where we have interchanged the raising and lowering of the indices on the left 
hand side. The equations can now be solved for gradients of f by inserting 
Eq. (|74D into ([T^): 



u^d^v'^ = vlvo,d,v'^. (7^ 
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Using relation (^) to replace the velocity field by the derivative of its po- 
tential $, and rearranging terms, we get 



d, - ) d>^<l> = 0, (79) 



which can alternately be written as 



((7% - = 0, (80) 



or 



where we denote the tensor ^ = -^v^Vq^ and ^ = ^ — ^ 



Before resuming with an interpretation of the result (|80|), we wish to 
express the variation in density of events n', in terms of the potential $. 
Inserting relation (^) into Eq. (^) and rearranging terms we get 

9^ [n' - = 0, (82) 

implying that the term in parenthesis is constant. By definition though, any 
constant part of n' is to be included in no, so the additive constant must be 
identically zero and we are left with the following relationship 



n' = ^v^^d,^. (83) 
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We now return to Eq. (pOf), and assume for simplicity that the background 
fluid is stationary i.e. Vq = {vq, 0, 0, 0) is pure time-like and constant. We 
then get 



d^d^ +lJ^\dod^ $ = 0, 




(84) 



or, on separating time and space derivatives, 



V' 



2 



( 




(85) 



This is a wave equation for which the modified sound velocity is 



cu 



(86) 



u = 



,02 



and we can investigate the various limits as follows. When we take the fluid to 
its noncompressible limit, or mathematically u —>■ oo, we get for the sound 
velocity — > c. This agrees with our conclusions of the previous section, 
namely that a noncompressible fluid gives rise to light speed sound waves. 

On the other hand if we take the fluid to be 'sufficiently' compressible, 
or in other words u ^ vq, the velocity becomes linear in u, u ^ -^^u. This 
is feasible for in the nonrelativistic limit when vq c, we then get the usual 
nonrelativistic result for compressible flow namely, u ^ u. These results are 
summarized in figure |I|. 



32 



5 Viscous fluids 

Until now we have dealt solely with "ideal fluids" which do not sustain shear 
stresses during motion. The justification for our previous equations was 
that for many materials the shear stresses occurring during motion are small 
compared with the pressure. In practice though, all liquids and gases are in 
fact able to sustain shear forces, and we now wish to take these into account 
in a consistent covariant manner. In order to obtain equations describing the 
motion of a viscous fluid, we have to include some additional terms in the 
constitutive equations of the fluid. We therefore write the stress tensor a^" 
in the form 

a^""" = -pg^""" + a^"". (87) 

The second term is the extra stress tensor or viscosity stress tensor resulting 
from frictional forces between different layers of fluid. We can, with some 
restrictive assumptions, establish a general form for the tensor . If we 
assume that our fluid does not "remember" its past history or initial con- 
flguration but rather that its motion is governed solely by the immediately 
preceding state, then ct'^'^, must depend only on the velocity held of the fluid. 
Processes of internal friction occur in a fluid only when different fluid parti- 
cles move with different velocities, so that there is a relative motion between 
various parts of the fluid. On the other hand if the distances between ad- 
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jacent parts of the fluid are kept constant during the motion and assuming 
that forces between events are a function of the distance separating them 
Ax^Ax'^, then there will be no friction. In other words a frictional dissi- 
pation could arise only in those regions of the fluid continuum undergoing 
distortion, excluding places at which the fluid moves uniformly as a rigid 
body. Hence a^'^ must depend on derivatives of the velocity field rather than 
the velocity field itself. When the velocity gradients are small we can, to 
some approximation, suppose a^'' to be a linear function of the derivatives 
d^v'^ omitting higher orders of powers and derivatives. There can be no terms 
independent of d'^v'^, since a^'^ must vanish for constant v'^. 

A rigid rotation of the fluid in space time must also be excluded from 
affecting the viscous force, and we shall show that this implies that only the 
symmetric combination of derivatives d'^v'^ + d'^v'^, can be contained in a'^". 
We must note here that by rigid rotation in spacetime, we mean the group 
of transformations A'^'^, which keep the covariant distance dx^j^dx^ constant. 
This is precisely the group of pure Lorentz transformations, which fulfill the 
requirement 

A^^A^^'^g^^. (88) 
We can show quite generally that for an infinitesimal transformation A*^ i, — 
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y + e^ jy, the infinitesimal displacement tensor y, must be antisymmetric, 
as follows: 

= g^'yg^-' + e^yg,'' + g^ye^'' + 0{6^) 

^ g% + {e% + e,''), (89) 

where in the third line we have neglected terms of second order in e. Com- 
parison of Eq. (pop with ( pSD yields 

e^'y = -e^>', (90) 

proving the antisymmetry of e'^^. 

Consider now a small region in the flow around a point x^, at which the 
fluid can locally be considered to be rotating rigidly (figure 0). Our claim 
is that at there are no frictional forces within the fluid and therefore the 
the viscosity stress tensor must vanish at that point. Let y represent the 
angular velocity of the fluid rotating in the /i — z/ plane. The velocity field 
round will then by definition be 

V^=^^yX\ (91) 
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where x'^ here and in the following discussion is the displacement from the 
center of rotation x^. Now consider a fluid particle at some place being 
rotated by an infinitesimal amount to a new position x''^. There are two 
relations connecting x and x': 



x'^ = A'^^x^ (92) 
x"' = x^ + Stv", (93) 



where 



A^. = g^. + e''., (94) 

is an infinitesimal Lorenz transformation. Inserting Eq. (^) into (^3]), and 
comparing with Eq. ( ^2]) results in 

A'', = g% + STn% (95) 



Using relation (Q) we can get a form for 17'^ ^ in terms of ^ 



brQ}', = e^,, (96) 



proving that Vt^ ^ too, is antisymmetric. Taking derivatives of relation (91) 



while considering the fact that fi'^ y is antisymmetric shows that the combi- 
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nation d^v'^ + d'^v^ vanishes for a rigidly rotating fluid whereas the antisym- 
metric combination d^v'^ — d'^v^, does not. As the viscosity stress tensor a^'^ 
must vanish for rigidly rotating motion, we deduce that it must contain just 
the symmetric combination of derivatives. 

The most general tensor of rank two satisfying all the above conditions is 

a>"' = aidf^v" + d'v^") + hg^'^dy, (97) 

where a and h are independent of the velocity. It is convenient, however, to 
replace a and h by other constants and write the equation in another form 
which lends itself more readily to interpretation 

= ^(^My^ + _ ^gi^^Q^yi) + Cg^^'dy. (98) 

We shall call the constants 7] and C, the shear viscosity and bulk viscosity 
respectively. The terminology derives from considering the two following 
basic flows (see figure |^): 

1. For a "pure" shear flow = u{n'y)n'^, where and n''^ are mutually 
orthogonal unit vectors so that the direction of the velocity gradient n'^ 
is perpendicular to the direction of the velocity, it can be shown that 
only the first term in a^" survives, leaving o"^^ = ur]{n'^n"^ + n'^n'^). 
So the coefficient t] alone governs the resistance of the fluid to shear 
distortions. 
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2. For a "pure" radial flow = (m/4)x^ in which the material expands 
in a radial direction, the the first term in Eq. ( P5| ) vanishes leaving 
^fiu _ y^Qg^^^ ^ which shows that C is proportional to the resistance of 
the fluid to expansion of bulk. This can be understood also by noticing 
that the expression in parentheses in Eq. (|98|) has the property of 
vanishing on contraction with respect to /i and v. 

The equations of motion can now be obtained by adding the expressions 
dyCF^'' to the left hand side of the nonviscous flow equation (^). Thus we 
have 



(99) 

This is the most general form of the equations of motion of a viscous fluid. 
If we assume, however, that the viscosity coefficients do not change notice- 
ably throughout the fluid, then they may be regarded as constant and can 
therefore be taken outside the gradient operators. We then have 



- d^p + r]d,d''v^' + (C + \v) d^dy = n + v'dy^ . (100) 

Further simplification can be obtained by assuming the bulk viscosity ( small 
compared to the shear viscosity rj. 
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If the fluid is incompressible then di,v 
hand side of Eq. ( |10UD vanishes giving 



= and the third term on the left 



- df'p + r]dud''v'' = n + v^'dyV^j . (101) 

On solving a specific dynamical problem one must also write down the 
boundary conditions for the equations of motion of a viscous fluid. We as- 
sume that there are friction forces acting between the boundary surface and 
the fluid such that the layer of fluid immediately adjacent to the boundary 
is brought to complete rest. Accordingly, the boundary conditions on the 
equations of motion require that the fluid velocity should vanish at fixed 
solid surfaces: 

= 0. (102) 

We emphasize that for a viscous fluid both the normal and tangential velocity 
components must vanish, in contradistinction to ideal fluids for which it is 
required only that the normal component vanish. In a general case with 
boundaries moving in r, the velocity of the fluid at the boundary should be 
equal to the velocity of the moving surface. 
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6 Conclusions 

We have studied a continuum flow of events in space-time parameterized 
by r, tlie invariant, universal, liistorical time. In tfie process of evolution, 
as r changes uniformly, the events move in space-time, generating a dense 
continuum of world lines which constitute a physical flow of particles. The 
velocity field is in general a function of both space-time and r and is governed 
by equations of motion and continuity. The solution of the equations for the 
velocity field provides, in principal, a solution to the physical problem, but 
is in general rather complicated. 

We considered an inviscid fluid and derived an analogue of the Bernoulli 
equation. In the special case of irrotational flow we proved that the velocity 
field can be derived by a potential. If in addition, the fluid is incompressible, 
it has been shown that the potential must obey the d'Alembert equation, 
and thus the problem is reduced to solving the d'Alembert equation with 
specific boundary conditions. The solutions consist of a background flow 
(The nonrelativistic solution) over which wave like ripples propagate with 
group and phase velocities equal to that of light. This is to be compared 
with Synge , who actually defines noncompressibility in such a way as to 



achieve this result. 

Relaxing the constraints and allowing the fluid to become linearly com- 
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pressible one can derive a wave equation, from which the sound velocity can 
again be computed. For a stationary background flow, it has been demon- 
strated that the sound velocity, attains its correct values for the incompress- 
ible and nonrelativistic limits. It is in principal possible, under some re- 
strictions, to compute the sound velocities also for a non stationary and 
non uniform background flow. Precise computations and predictions in this 
direction are yet to be studied. 

Finally viscosity was introduced and the use of some general arguments of 
symmetry and isotropy enabled the formulation of equations for the motion 
of a viscous fluid. Only introductory comments were made about this topic 
and the equations derived are yet to be applied to specific problems. 

It must be emphasized that the velocity and density (and any other) fields 
that we mention, describe quantities related to the flow of events at a par- 
ticular time T. In order to derive a physically observable quantity related to 
particles, an averaging must be performed over r. The method we use is that 
of integrating over r with the density of events as a weight function. An ex- 
act transcription from event properties to world-line, or particle, properties, 
however, is not always possible due to nonlinearities of these quantities. Al- 
ternatively the observed quantities can be extracted from the zero frequency 
(in r) components of the fields. 
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Figure 1: a; as a function of the inverse compressibility u, and energy vq 
respectively. 
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Figure 2: Infinitesimal rotation of fluid. 
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Figure 3: Two basic distortions: left: Radial flow, right Shear flow. 
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a; as a function of the inverse compressibility u, and energy vq respectively. 



Infinitesimal rotation of fluid. 



Two basic distortions: left: Radial flow, right Shear flow. 
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